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^  Abstract 

The  evolution  of  unsteady  boundary  layers  on  the  line  of  symmetry  of  a 
thin  prolate  spheroid  In  uniform  motion  at  constant  angle  of  attack  after  an 
Impulsive  start  Is  studied  for  a  prescribed  pressure  distribution  and  results 
have  been  obtained  for  angles  of  attack  ranging  from  30*^to  50^7  The  study 


shows  that  if  ot  <  acM!°)  and  the  steady-state  boundary  layer  exists  every¬ 
where  on  this  line  of  symmetry,  the  unsteady  solution  approaches  Itk  steady- 
state  solution  rapidly  without  any  significant  special  features,  -^he  same  is 
true  for  a  <  a  on  the  windward  side:  on  the  leeward  side,  howevef,  the 

c  V 

steady-state  boundary  layer  separates  and  the  solution  terminates.  ^  The 
unsteady  boundary  layer  is  Initially  unseparate d  but  develops  a  region  of 
reversed  flow  after  a  finite  time.  A  short  time  later,  the  streamwise  dis¬ 
placement  thickness  develops  a  pronounced  peak  and  leads  to  a  singularity  of 
the  same  type  as  that  observed  by  van  Dommelen  and  Shen  on  a  circular  cylinder 
started  impulsively  from  rest. 


1.0  INTRODUCTION 

Our  knowledge  of  unsteady  boundary  layers  has  been  enhanced  in  recent 
^  "years  by  a  combination  of  experimental  and  computational  investigations.  The 
L  former  have  been  concerned  mainly  with  the  problem  of  airfoils  which  oscillate 
L„.1n  subsonic  flows  but  include  a  small  number  of  transonic  flow  investigations. 
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The  latter  have  Involved  the  solution  of  boundary-layer  equations  for  flows 
over  cylinders  and  airfoils  and  have  been  concerned  largely  with  the  nature  of 
the  singularity  which  occurs  near  separation.  Solutions  of  the  Navler-Stokes 
equations  have  also  been  obtained  to  provide  an  overall  understanding  of  the 
flow  patterns  at  conditions  near  or  at  dynamic  stall.  Useful  reviews  of  prev¬ 
ious  work  have  been  provided  by  Tellonls  (1979),  McCroskey  (1982),  Williams 
(1977)  and  Cebeci  (1982). 

The  experimental  Investigations  have  been  carried  out  mainly  to  improve 
understanding  of  flow  around  helicopter  rotor  blades  and  Include  the  low-speed 
measurements  of  McCroskey  et  al.  (1982),  Carr  et  al.  (1977),  Carr  and  McAlister 
(1983),  Young  (1982),  Geissler  (1983)  and  Cousteix  et  al.  (1981)  and  transonic 
flow  measurements  of  Davis  and  Malcolm  (1980)  and  Dowell  et  al.  (1981).  The 
range  of  measurements  of  McCroskey  et  al.  and  Carr  et  al.,  in  particular,  is 
extensive  and  includes  more  than  fifty  combinations  of  Mach  number  and  param¬ 
eters  of  the  unsteady  motion  for  each  of  eight  airfoil  sections.  As  a  conse¬ 
quence,  four  flow  regimes  have  been  identified  and  correspond  to  ho  stall, 
stall  onset,  light  stall  and  deep  stall.  It  appears  that  the  breakdown  of  the 
unsteady  boundary  layer  leads  to  a  large  vortex  which  Is  formed  near  the 
surface  at  large  angles  of  attack  and  causes  stall  to  occur  shortly  thereafter. 

Computational  investigations  of  oscillating  airfoil  flows  have  been  of  two 
main  types.  In  the  first,  solutions  of  the  Navier-Stokes  equations  have  been 
obtained,  for  example  by  Mehta  (1977)  for  incompressible  laminar  flows  and  by 
Shamroth  (1981)  for  compressible  turbulent  flows.  Further  studies  of  this 
type  are  clearly  necessary  and  results  so  far  show  qualitative  features  of  the 


flow  field  and  lift  curves  which  are  again  in  qualitative  agreement  with 
experiment.  The  second  approach  has  Involved  the  solution  of  the  boundary- 
layer  equations.  The  Investigations  of  Cebeci  and  Carr  (1981,  1983)  have 
reported  results  for  an  external  velocity  distribution  typical  of  those  found 


near  the  leading  edge  of  thin  airfoils.  The  existence  of  a  singularity  In  the  v  coj»s 
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solutions  Is  evident  In  the  vicinity  of  the  leading  edge  at  the  higher  angles 
of  attack  and  has  led  to  the  subsequent  Investigations  of  Its  nature,  see  for 
example,  Cebecl,  et  al.  (1983). 

The  need  for  more  fundamental  Investigations  of  the  use  of  boundary-layer 
equations  to  represent  unsteady  flows  Is  clear  from  the  airfoil  Investigations 
of  the  previous  paragraph.  In  this  connection,  the  oscillating  flat  plate 
Investigations  of  Coustelx  et  al.  (1981)  and  the  contributions  of  Telionis 
(1974),  van  Dommelen  and  Shen  (1982),  Smith  (1982)  and  Cebecl  (1979,  1982)  for 
flow  over  a  cylinder  Impulsively  started  from  rest  are  particularly  relevant. 

In  addition,  Williams  (1982)  and  Williams  and  Stewartson  (1983)  have  made 
Important  contributions  to  our  understanding-  of  the  nature  of  the  singularity 
and  its  consequences.  Perhaps  the  most  Inportant  contribution  has  been  that 
of  van  Donmelen  and  Shen  (1982a)  who  solved  the  boundary- layer  equations  In 
Lagranglan  form  and  revealed  the  existence  of  the  singularity.  This  result 
suggests  the  need  for  Interaction  between  the  viscous  and  Invlscld  flow 
equations. 

Though  a  number  of  studies  have  been  conducted  to  improve  our  understand¬ 
ing  of  unsteady  two-dimensional  boundary  layers,  very  little  work  has  been  done 
for  unsteady  three-dimensional  flows.  Experimental  information  is  lacking,  but 
the  similarity  between  the  two-  and  three-dimensional  equations  suggests  that 
the  same  phenomena  may  occur.  The  three-dimensional  boundary  layer  on  a  body 
of  revolution  at  angle  of  attack  is  obviously  more  complicated  than  for  air¬ 
foils  or  cylinders  but  Is  similar  on  the  line  of  symmetry,  provided  that  cross- 
flow  gradients  are  taken  Into  account.  In  the  case  of  steady  flow  past  thin 
spheroids.  It  Is  known  [see  Cebecl  et  al.  1981a]  that  separation  does  not 
occur  on  the  leeward  side  until  well  past  the  maximum  thickness  if  the  angle  of 
attack  a  Is  less  than  a  critical  value  a„(«41°),  but  for  a  >  a  there  Is  a 
dramatic  change  with  separation  occurring  very  close  to  the  nose.  If  the 
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external  velocity  Is  prescribed,  this  abruptly  terminates  the  Integration.  For 
thin  airfoils  there  Is  a  parallel  situation  but  the  critical  angle  Is  now  of 
the  order  of  magnitude  of  the  airfoil  thickness. 

In  this  paper  we  report  a  study  of  the  unsteady  boundary  layer  on  the  line 
of  synmetry  to  determine  the  relation  between  unsteady  separation  and  singu¬ 
larities  In  the  solution.  For  the  first  step  we  concentrate  on  the  separation 
problem  and  do  not  consider  the  effect  of  moving  stagnation  points  which  Is 
Important  but  adds  considerably  to  the  computational  difficulties.  The  par¬ 
ticular  problem  we  study  Is  the  development  of  a  boundary  layer  on  a  thin 
spheroid  In  uniform  motion  at  constant  angle  of  attack  after  an  impulsive 
start.  Angles  of  attack  ranging  from  30°  to  50°  are  chosen  for  this  purpose. 

Of  these  a  *  45°  corresponds  to  nose  separation  which  Is  marginal  while  a 
=  50°  corresponds  to  a  strong  steady-state  singularity.  The  problem  has  been 
formulated  for  the  general  case  of  spheroids  at  angles  of  attack  and  special¬ 
ized  to  thin  bodies,  and  specifically  those  for  which  the  thickness  ratio  t 
Is  vanishingly  small. 


2.0  APPROACH 

With  respect  to  a  cylindrical -polar  set  of  coordinates  (x,r,9)  and 
origin  0,  let  the  equation  of  the  spheroid  be 

A2  +  r2  =  t2  (1) 

where  t  Is  a  positive  number  less  than  unity,  which  may  be  regarded  as  a 
measure  of  the  thickness  ratio  of  the  prolate  spheroid.  Here  we  have  nondi- 
menslonallzed  the  coordinates  and  time  t  taking  the  semi-major  axis  and  the 
velocity  at  Infinity  to  be  unity.  We  assume  that  at  an  Infinite  distance  from 
the  spheroid  the  velocity  of  the  fluid  is  In  a  direction  lying  in  the  merid¬ 
ional  plane  9  =  0  and  making  an  angle  a  with  Its  major  axis.  Further  we 
assume  that  outside  the  boundary  layer  the  invlscid  flow  Is  Irrotatlonal  and 
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attached  and  neglect  the  circulation  around  the  spheroid.  Although  the  circu¬ 
lation  Is  generally  Important,  particularly  at  large  angles  of  attack,  near 
the  nose  Its  effect  Is  mainly  to  change  the  value  of  a  and  so  the  general 
conclusions  of  this  Investigation  may  be  modified  to  take  account  of  It  In 
defining  a  In  terms  of  the  position  of  the  forward  stagnation  point.  The 
achievement  of  our  goal  of  determining  the  development  of  the  boundary  layer, 
given  the  main-stream  velocity,  would  be  an  important  step  forward  In  the  task 
of  predicting  the  circulation  around  bodies  of  revolution  at  Incidence. 

With  these  assumptions  the  velocity  of  slip  on  the  spheroid,  according  to 
Invlscld  theory,  has  components 

ue  3  V„(t)  coso  cosp  +  Vqn(T)  sina  sin0  cos9 

(2) 

we  *  Vgg(T)  Sina  Sln0 

where 


V0(T) 


(1  -  t2)3/2 


n  .  Ty/2  -  1/2  x2log[(l  +  /l  -  t2)/(1  -  /l  -  T2)] 
ZV0(tJ 

V90(t)  *  ?V0Tt)  ~T 


(3) 


cosB  = 


(1  - /)T/2 
[1  +  XZ(T2  -  1)]1/2 


(4) 


and  a  Is  a  specified  function  of  time,  being  constant  In  the  present  study. 
The  configuration  in  the  meridional  plane  is  Illustrated  In  figure  1  with  S 
the  stagnation  point  of  the  invlscld  flow  defined  by  0  =  0,  x 
■  -{1  -  t2  tan2a} T//2. 


Turning  to  the  boundary  layer  we  define  (u,v/v,w)  to  be  the  velocity 
components  respectively  along  the  meridional  lines  0  *  const,  along  the  nor¬ 
mals  to  the  spheroid  and  In  the  azimuthal  direction  where  v  Is  the  kinematic 
viscosity  of  the  fluid  and  Is  assumed  small.  Then 
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(5) 


fx  *h2u*  +ls'  *hlw*  +!y  *h1h2v*  a  0 

“  2 

iU  .  U  iU  ,  W  iU  .  JU  .  A  ,  1  3P  .  3  U 

Tf  +W7?x  *iq^5+  v3y  +  ™2  *  "  TTfSx l6) 

»  y 

3u  ,  u  3w  w  3w  3w  ~  !_lR  +  iL£  m 

H  +  K737  +  H71e  +  V1y  uw,!2  HJ3t  +  ^  m 


where  jvv  measures  distance  along  the  outward-drawn  normal  from  the  body. 


1  +  x2(t2  -  1)  1/2 

h,  -  t-  -  ■-%  -’l  . 

1  .  1  -  xz 

are  the  metric  coefficients, 


h2  .  t(  1  -  x2>1/2 


K2  ■[i;7(7"-ljfrr~x^ 


(8) 


(9) 


is  the  geodesic  curvature  of  the  surface  lines  x  *  const,  and  p  Is  the  reduced 
pressure.  The  appropriate  boundary  conditions  are 


u  *  v  =  w  *  0  at  y  *  0 


(10a) 

u  -*■  ue,  w  ♦  we  as  y  +  a> 

The  specification  of  the  problem  Is  completed  by  assigning  initial  conditions 
on  u,w  with  respect  to  space  and  time.  The  spatial  conditions  are  that  on  the 
line  of  symmetry  w  »  0  and  u  vanishes  for  all  y  at  that  value  of  x  for  which 
ue  »  0.  This  Is  physically  and  mathematically  sensible  for  impulsive  prob¬ 
lems  in  which  the  stagnation  point  Is  fixed  but  Is  not  necessarily  correct 
when  the  angle  of  attack  Is  varying  with  time,  see  Cebeci  and  Carr  (1981)  for 
example.  The  temporal  conditions  are  that  at  t  *  0  both  u  and  w  vanish  on  the 
body  y  a  0  on  the  line  of  symmetry  but  are  equal  to  their  external  values 
elsewhere.  Thus 

u  ■  w  *  0  at  y  *  0,  t  *  0  for  all  x 


while 


(10b) 


u  *  u„,  w  »  w„  for  y  >  0,  t  *  0  and  all  x. 
e  e  v 
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In  the  early  studies  of  the  properties  of  the  boundary  layers  on  prolate 
spheroids  conducted  by  Wang  (1970)  and  Hirsh  and  Cebecl  (1977),  some  diffi¬ 
culty  was  experienced  In  continuing  the  solution  past  the  nose  at  x  ■  -1,  pri¬ 
marily  because  of  the  singularities  In  the  properties  of  hj,  h2  and  Kg  there. 

A  common  procedure  was  first  to  perform  the  Integration  along  the  line  of  sym¬ 
metry  from  the  stagnation  point  to  as  near  the  nose  as  possible  and  then  to 
jump  across  the  nose  to  the  same  value  of  x  on  the  leeslde  (0  *  n)  assuming 
that  the  flow  properties  are  essentially  unchanged.  After  that  Integration  on 
the  leeslde  may  be  continued  as  far  as  separation.  Afterwards  the  procedure 
may  be  extended  to  more  general  points  In  the  neighborhood  of  the  nose.  This 
technique  is  effective  at  moderate  values  of  t  [t  *  1  corresponds  to  a  sphere] 
but  leads  to  difficulties  as  t  0  especially  at  high  angles  of  attack.  Cebecl 
et  al.  (1980)  demonstrated  that  the  singularity  may  be  removed  by  a  suitable 
transformation  of  the  surface  coordinate  system  enabling  a  smooth  passage  to 
be  made  around  the  nose.  The  transformation  Is  equally  effective  for  unsteady 
boundary  layers  and  Is  now  explained  for  the  limiting  case  of  a  paraboloid 
which  corresponds  to  t  *  0. 

The  first  step  Is  to  define  new  surface  coordinates  by 
dS  _  hldx  _  [1  +  x2(xZ  -  1)]1/2  . 

s  -r -  a  - 3 -  ClX 

5  n2  t( 1  -  X2) 

(ID 

X  «  S  cose,  Z  *  S  s1n0,  Y  *  y/r 
with  S  »  0  at  x  *  -1 ,  and  new  velocity  components  by 

u  *  U  cos0  +  W  sin0,  w  *  W  cos9  -  U  sin©,  v  ■  V/t  (12) 

The  purpose  of  this  transformation  Is  to  convert  the  polar  form  of  the  equa¬ 
tions  (5)  -  (7)  Into  a  quasl-rectangular-Carteslan  form  which  Is  free  of 
singularities.  The  governing  equations  reduce  to 

N(|j|  +!l)  +!t  "  L(UX  +  WZ)  *  0  {13) 

7 


t 


(14) 


||  +  N(u  |U  +  w  |U)  +  LW(WX  .  uz)  +  V  U  •  &,  + 


32U 

a? 

32W 


||  +  N(U  ||  +  W  ||)  -  LU(WX  -  UZ)  +  V  ||  »  %2  +  ^jjf  (15) 

where  %2  are  Pressure-gradient  parameters  defined  by  the  requirement 


that  u  -►  ua,  w  w.  as  Y  +  «  and 
e  e 


St 


ITT 


L  s  r  +  k2) 


For  the  paraboloid  we  let  t  -*■  0  after  defining 

n,ihM 


(16) 


(17) 


when 


S  . - R— - exp  [(1  +  p2)1/2  -  13 

(1  +  p*)vz  +  1 


"■f  * 


,  (1  +  P2)T/2  -  1 

pS(l  +  p2)172 


(18) 


(19) 


and  so  Eqs.  (13)  -  (15)  are  independent  of  t.  We  write  (Ue>  Wg)  for  the 
limits  of  (ue  cose  -  wfi  sine,  ug  sine  +  wg  cos  e)  as  t  -*•  0 


and  then 


so  that 


U I  =  Px  cosa  _  2(1  _  sina 

e  S(1  +  p2)172  S 


W  -  PZ  ■c-g5T7y  +  s1 

e  S(l+p2)1/Z  S 


not 


.  3U  3U 

Jt,  -  «(U,  ^  ♦  We  ♦  U.IW.X  -  UeZ) 


(20) 

(21) 


(22) 


3W  3W 

h  ■  N<“e  7T  *  -  LW  -  V> 


(23) 


These  equations  are  free  of  all  geometric  singularities  and  In  particular 
at  the  nose  p  ■  0.  There  Is  therefore  now  no  special  problems  about  Integrat¬ 
ing  the  equations  through  the  nose  although  It  should  be  noted  that  the  equa¬ 
tions  are  only  appropriate  at  distances  from  It  cor  espondlng  to  p  ■  0(1). 

This  Is  the  natural  length  scale  for  the  paraboloid  but  In  terms  of  general 
thin  axlsymmetrlc  bodies  this  distance  Is  0(t  )  from  the  nose. 


2.1  Llne-of -Symmetry  Equations 

These  may  be  deduced  from  the  general  form  Eq.  (13)  et  seq.  given  in  the 
previous  section  using  similar  arguments  to  those  In  Cebeci  et  al.  (1981).  We 
shall  not  repeat  them  In  detail  but  merely  state  the  transformations  and  the 
final  form.  First  we  write 

U  »  Uq(X,Y,T)  +  0(Z2),  V  *  V0(X,Y,T)  +  0(Z2)  (24) 

W  -  Z  exp[l  -  /l  +  p2]Wj(X,Y,T)  +  0(Z3)  (25) 


and  allow  for  negative  values  of  X  by  permitting  p  to  become  negative.  When 
p  <  0  the  sign  of  S  in  Eq.  (18)  must  be  changed  and  generally  X  *  S  sgnp  in 
the  limit  Z  0.  On  the  leeside  line  of  symmetry  p  <  0. 

We  now  substitute  Eqs.  (24),  (25)  Into  (13)  et  seq  and  take  the  limit 
Z  0.  The  external  velocity  components  reduce  to 

2  1/2 

„  _  p  cosa  -  2  slna  ,,  C(1  +  p  )  +  l]cosa  +  2p  sina  ,„,.N 

«  (7.7)1'*  ’  16 - o  ♦  ’ 

To  put  the  equations  into  the  most  convenient  form  for  numerical  integra¬ 
tion  of  the  impulsive  problem,  with  primes  denoting  differentiation  with 
respect  to  n  we  write 


- ,  T  -  *  p2),/4v0  -  ”, 

(l+p2)T/V/2  0  0  2(1  +  p2)  ' 


(27) 
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where  f  and  g  are  functions  of  p,n»t.  Then 


dU_ 


Yo  "  *  Uoe  Ip  '  f(a1Uoe  +  “3T  )  "  *29 


where 


(i  ±  p2)1/2 


1  2d  +  pz)  (pz  +  d1^  + 1  *  2  i  +  (T77?77 

and  the  appropriate  boundary  conditions  for  f  and  g  are: 

f  •  f1  s  g  8  g'  a  0  at  n  ■  0, 


(29) 


(30) 


(31) 


f  ♦  1,  g'  w 


le 


as 


The  momentum  equations  reduce  to 


f-  -  p,ff  ♦  P6gf  ♦  p2[l  -  (f)2]  *  p,3rf-  -  r,  If  ♦  P)0(f  If  -  f’  If 

(32) 

g-'  +  P-jfg"  +  P6gg"  +  P3(f')2  +  P4(g’)2  +  Pgf'g'  +  P12  +  P13ng" 
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(33) 
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2,1/2 


(34) 


These  two  equations  (32)  and  (33)  are  In  a  convenient  form  for  numerical 
Integration  since  at  t  =  0  f,g  and  their  derivatives  are  given  by 


where 


2  2 

f  *  nerf(Xn)  — — (1  -  iX  n  ), 

Xy'-rr 
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«  ■  w„f 

g'  ■  Wlef'  (35) 

-  WIef 


Further,  the  Integration  In  the  direction  p  Increasing  may  be  assumed  to  start 
at  the  stagnation  point  p  *  2  tana,  where  f,g  are  Independent  of  p  and  sat¬ 
isfy  partlal-dlffer* ntial  equations  in  n,t  only,  and  proceeds  in  a  straight¬ 
forward  way  since  Uoe  >  0  when  p  >  2  tana.  A  similar  remark  applies  to  the 
integration  in  the  direction  of  p  decreasing  since  UQe  <  0  when  p  <  2  tana. 


3.0  SOLUTION  PROCEDURE 

The  numerical  solution  of  equations  (32)  and  (33)  was  obtained  using  the 
Box  method,  a  two-point  finite-difference  method  developed  by  H.  B.  Keller  and 
extensively  used  for  two-dimensional  flows,  both  steady  and  unsteady,  and  for 
three-dimensional  steady  flows.  Here  we  use  two  versions  of  the  Box  method, 
the  regular  box  in  regions  where  f 1  ^  0  across  the  layer  and  the  zig-zag  box 
In  regions  where  f 1  <  0  at  any  value  of  n.  Specifically  in  advancing  the 
solution  to  a  new  value  of  t  or  p,  this  property  of  f '  determined  the  choice 
of  the  box  scheme.  The  details  of  the  box  scheme  are  described  in  Bradshaw  et 
al.  (1981). 
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Calculations  of  the  unsteady  boundary  layers  were  carried  out  for  angles 
of  attack  a  of  30°,  40°,  45°,  50°.  Provided  separation  did  not  occur,  the 
solution  of  the  equations  by  the  regular  box  was  appropriate,  no  numerical 
difficulties  were  encountered  and  the  solution  was  smooth.  The  step-sizes 
chosen  for  all  computations  were  Ap  »  0.2,  At  ■  0.2;  In  selected  cases  compar¬ 
isons  were  made  with  studies  using  step-sizes  of  Ap  *  0.1  and  At  ■  0.1  or 
smaller  and  the  differences  were  negligible.  However,  once  separation  occur¬ 
red,  the  regular  box  became  prone  to  Instabilities  and  was  replaced  in  this 
domain  by  the  zig-zag  box  over  all  n  and  the  step-sizes  In  p  and  t  were 
reduced.  Various  calculations  were  carried  out  with  a  nonuniform  mesh  In  p,  t 
In  the  neighborhood  of  separation.  The  smallest  value  of  Ap  was  0.0005  near 
the  separation  point;  the  step  size  in  t  was  progressively  reduced  as  the 
calculations  proceeded,  the  smallest  value  being  0.01.  Elsewhere  Ap  was 
fixed  at  0.2.  Comparisons  with  solutions  using  larger  step-sizes  gave 
confidence  that  the  results  were  reliable  outside  the  separated  region  for  t 
<  7,  but  that  inside  the  separated  region  smaller  step-sizes  In  p  and  t  were 
necessary  to  avoid  small  oscillations  which  developed  for  t  >  5.5  at  a  * 

50°. 

4.0  RESULTS 

The  principal  properties  of  the  unsteady  boundary  layers  are  the  skin- 
friction  components  defined  by 

3U 

Ug(o,p,t)  =  t"1/2Ue  f»(o,p,t)  =  (1  +  p2)1/4  (j^)w  (36a) 

W^(o,p,t)  «  f1/2g«(o,p,t)  =  (T  +  p2)T/4  (-^r)w  (36b) 

and  the  displacement  thicknesses  defined  by 

12 


Aj(p,t)  -  Ue  t'fC  11m  (n  -  f)  -  (1  +  p2)"1/4  /  (U  -  U  )dy 


(37a) 


0  X  ® 


A2(p,t)  .  tT/2  11m  [nWg  -  g)  -  (1  +  p2)“1/4  J  (Wel  -  W,)dy  (37b) 

These  are  displayed  In  Fig.  2  for  a  ■  30°  when  separation  does  not  occur. 

In  Fig.  3  we  display  the  properties  of  UQ  for  a  *  40°,  a  solution  which 
Is  on  the  verge  of  separation,  and  In  Fig.  4  we  show  all  the  functions  [Eqs. 

(36)  and  37]  for  a  *  50°  when  separation  clearly  occurs.  Where  appropriate 
the  corresponding  steady-state  results  of  Cebecl  et  al.  (1980)  are  included  In 
these  graphs  for  comparison  purposes. 

The  graphs  show  that  for  p  >  0,  l.e.  on  the  windward  line  of  symmetry, 
the  steady-state  solution  Is  rapidly  approached  as  t  Increases  and  is  essen¬ 
tially  established  for  t  >  1. 

For  p  <  0,  l.e.  on  the  leeward  line  of  symmetry,  the  flow  properties  In 
the  steady  state  have  two  Interesting  features.  First  for  any  value  of  a, 
the  cross-flow  displacement  thickness  &2  rapidly  decreases  with  p  and  soon 
becomes  negative,  even  If  separation  has  not  occurred,  and  In  fact  Is  exponen¬ 
tially  large  and  negative  when  p  Is  large  and  negative.  The  unsteady  data  for 
a  =  30°,  Fig.  2d,  confirms  this  trend  but  the  limits  t  ®,  p  ®  are  not 
Interchangeable,  i.e.  for  any  t  the  graph  of  flattens  out  for  large  enough  -p 
and  for  any  p  <  -2  the  graph  of  A,,  eventually  approaches  the  steady  state. 

I 

Second  the  steady-state  meridional  component  UQ  of  the  skin  friction  does 
not  tend  monotonlcally  to  a  limit  as  p  -*•  -  »,  but  even  at  small  values  of  a 
has  a  trough  followed  by  a  peak  due  to  the  overshoot  in  the  external  velocity  (26), 

Further,  at  a  >  ac  *  41°,  UQ  actually  vanishes  at  p  ■  pc(a),  when  p£ 

Is  a  point  near  p  *  -1  depending  on  the  particular  choice  of  a.  The  solution  is 
singular  at  this  point  and  the  calculation  must  terminate. 
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The  unsteady  solution  approaches  the  steady  state  as  t  ♦  ®  for  a  <  ac 
(Figs.  2  and  3)  but  more  slowly  on  the  windward  side.  Some  overshoot  may 
occur  but  It  Is  small  and  may  be  an  Indication  of  a  fall-off  In  accuracy  of 
either  the  steady-state  or  the  unsteady  solutions.  For  a  >  ac  (Fig.  4),  the 
same  remarks  apply  as  long  as  p  >  pc(a)  and  flow  reversal  has  not  occurred. 

For  smaller  values  of  p,  reversed  flow  occurs  and  there  Is  no  corresponding 
steady-state  solution.  If  o  *  45°,  separation  first  occurs  at  t  s  5.64  when 
p  *  -1.2;  It  remains  confined  within  the  range  -1.0  >  p  >  -1.4  and  Is  weak  until 

I 

t  »  6  when  the  computation  was  terminated.  It  might  be  argued  that  and 
^2  are  approaching  limit  states  as  t  -*■  «  but  Aj  Is  definitely  showing  signs  of 
a  pronounced  negative  minimum  near  p  *  -1.2. 

Separation  is  marginal  for  a  «*  45°  and  so  a  more  extensive  computation 
was  carried  out  for  a  »  50°.  In  this  case  separation  set  in  when  t  *  3.71 
at  p  *  -1.083  and  gradually  spread  out  as  t  increased  to  extend  around  the 
range  -0.908  <_  p  <_  -1.396  at  t  *  6  when  the  calculations  became  somewhat 
dubious  due  to  oscillations  and  Instabilities,  a  behavior  previously  observed 
In  relation  to  the  circular  cylinder  started  impulsively  from  rest,  Cebeci 

I 

(1982).  Over  the  same  range  of  t  the  variation  of  and  Wj  may  appear  to 

I 

remain  smooth  but  UQ  and  A  ^  develop  sharp  minima  just  downstream  from 
separation. 


5.0  THE  SINGULARITY? 

We  now  Investigate  whether  the  results  obtained  for  a  =  50°  allow  us  to 
conclude  that  when  the  computation  terminated  the  flow  properties  are  consist¬ 
ent  with  a  singularity  centered  just  downstream  of  separation  In  the  neighbor¬ 
hood  of  t  *  6. 

Two  possible  structures  for  a  terminating  singularity  to  the  solution 
have  been  proposed  of  which  one  Is  due  to  Smith  (1982).  In  his  model,  the 
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viscous  terms  are  significant  at  breakdown,  which  Is  a  generalization  of  the 
Goldstein  singularity,  and  both  skin-friction  and  displacement  thickness  are 
Infinite.  No  example  of  such  a  structure  has  yet  been  found  and  although  we 
cannot  rule  out  Its  relevance  In  the  present  flow,  our  Inclination  Is  to 
prefer  the  second  possibility,  due  to  van  Dommelen  and  Shen  (1982a),  In  which 
the  skin  friction  remains  finite  while  the  displacement  thickness  becomes 
Infinite.  A  full  description  of  this,  the  van  Dommelen  and  Shen  (1982b) 
singularity,  has  not  yet  been  given  and  some  questions  are  still  open  but 
strong  evidence  In  favor  of  Its  being  appropriate  to  an  unsteady  boundary 
layer  on  a  circular  cylinder  has  been  provided  by  van  Dommelen  and  Shen 
themselves  and  by  Cowley  (1983)  and  Ingham  (1983)  and  on  an  oscillating 
airfoil  by  Cebecl  et  al.  (1983). 

The  essence  of  the  singularity  is  Invlscid  and  independent  of  the 
pressure  gradient.  Thus,  in  the  case  of  a  two-dimensional  boundary  layer, 
with  equations 


3ji  av 
3x  ay 


0, 


(38a) 


au 

at 


+ 


2 

1  ap  .  a  u 

P  3x  +  ^2  * 


(38b) 


let  us  suppose  the  singularity  is  centered  at  t  *  t  ,  x  *  xQ.  Then  In  Its 
simplest  form  we  write 


u  *  -Uq  +  (tg  -  t)V2  u(x,y)  (39) 

when  0  <  t0  -  t  «  1,  |x  -  x0|  «  1, 

where 

x  ■  C(x  -  x0)  -  u0(t0  -  t)]/(t0  -  t)3/2,  y  »  y(t0  -  t)T/4 

and  the  determination  of  u  reduces  to  quadratures.  This  structure  must  be 
matched  to  other  structures  near  y  ■  0  and  as  y  -*•  ®  but  the  principal  gross 
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features  can  be  Inferred  from  (39).  The  displacement  thickness  has  a  peak  at 
^  -1/4 

x  »  0  which  becomes  Infinite  like  (tQ  -  t)  as  t  -*■  tQ.  On  the  other  hand, 
the  skin  friction  remains  finite  as  t  +  ts  for  all  x.  The  velocity  profile 
becomes  very  flat  as  t  ♦  tQ  provided  x  1,  with  u  ■  -uQ  over  a  length 
(tQ  -  t)  normal  to  the  surface,  the  adjustment  to  the  boundary  conditions  both 
on  the  surface  and  at  infinity  taking  place  over  effectively  finite  ranges  of 
values  of  y.  Further  details  of  the  structure  of  the  singularity,  on  the  basis 
of  a  Lagranglan  representation  may  be  found  In  van  Dommelen's  thesis  (1981) 
and  In  Elliot,  et  al.  (1983)  on  the  basis  of  an  Eulerian  representation  as 
exemplled  by  (38). 

This  theory  may  be  adapted  to  our  problem  at  least  in  principle.  We 
replace  u,x,y  In  (39)  by  f',p,n  and  neglect  the  left-hand  side  of  (32)  In 
the  neighborhood  of  p  *  p0,  t  =  t0.  Then  the  equation  for  t'  Is  the  same 
as  that  f or  u  and  we  may  Infer  that  f'  is  as  likely  to  develop  the  van 
Dommelen  and  Shen  singularity  as  Is  u  In  the  circular  cylinder  problem.  Thus 
the  displacement  thickness  A^  can  be  expected  to  develop  an  Increasingly 
sharp  peak  at  a  value  of  p  which  Is  a  linear  function  of  t  and  at  a  finite 
line,  a  needle-like  infinity.  Also,  the  corresponding  component  of  skin 

I 

friction  U  remains  finite  and  smooth. 

o 

The  evidence  in  Fig.  4  supports  these  conclusions.  A  more  severe  test  Is 
to  plot  the  streamwlse  displacement  velocity,  3A^/3p,  as  a  function  of  p  for 
various  t  and  this  is  done  In  Fig.  5a  after  smoothing  the  data  for  p  <  -0.95. 

It  Is  convenient  to  compare  these  graphs  with  one  (Fig.  6)  kindly  provided  by 
Cowley  for  solution  of  (38)  with  uQ  =  sinx.  In  which  the  singularity  occurs  at 
t  »  3.0.  The  similarities  in  form  are  striking;  the  principal  difference  being 
a  slight  drift  in  Mpn  3Aj/3p  away  from  separation  as  t  Increases  in  contrast 
to  the  prediction  in  (39)  [uQ  >  0].  The  value  of  p  Aj  drifts  towards  the 

I 

separation  point.  Further  the  more  detailed  plots  of  UQ  In  Fig.  5b  confirm 
the  smooth  behavior  of  this  function. 
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If  now  we  examine  (33)  we  see  that  the  same  form  for  g'  may  be  assumed 

•  i 

provided  gQ  Is  taken  equal  to  f  and  again  the  left-hand  side  of  (33)  may  be 

neglected.  This  means  that  A2  should  also  develop  a  needle-like  singularity  and 

Fig.  4d  shows  no  evidence  of  this.  However,  plots  of  3A2/3p  as  functions  of  p 

for  various  t,  displayed  In  Fig.  5c,  show  clear  signs  of  the  development  of  a 

strong  minimum,  similar  to  Fig.  6,  which  does  drift  towards  separation.  The 

alternative  possibility  that  A2  Is  finite  seems  to  require  that  as  n 

Increases  g‘  has  reached  the  value  W,  before  the  invlscld  zone  n  *  1  Is 

•e 

reached.  Were  this  the  case,  the  two  sides  of  (33)  would  not  balance  over  this  zone 

I 

since  then  f  “  fQ  i  1. 

6.0  DISCUSSION 

In  this  paper  we  have  examined  the  evolution  of  the  unsteady  boundary 
layer  on  the  line  of  symmetry  of  a  paraboloid  which  is  set.  Impulsively,  Into 
motion  at  t  *  0  with  uniform  velocity  and  at  an  angle  a  of  attack.  We  have 
shown  that  If  a  <  ac  (*41°)  and  the  steady-state  boundary  layer  exists 
everywhere  on  this  line,  the  unsteady  solution  approaches  It  reasonably  quickly 
and  without  any  significant  special  features.  The  same  Is  true  for  a  >  ac 
on  the  windward  side.  On  the  leeward  side,  the  steady-state  boundary  layer 
separates  and  the  solution  must  terminate  there  because  the  external  pressure 
gradient  Is  fixed.  The  unsteady  boundary  layer  Is  Initially  unseparated  but 
develops  a  region  of  reversed  flow  after  a  finite  time.  A  short  time  later  the 
first  displacement  thickness  A^  develops  a  pronounced  peak  and  we  advance 
arguments  for  believing  that  this  is  associated  with  an  Incipient  singularity 
which  terminates  the  calculation. 

The  results  obtained  here  are  also  relevant  to  other  bodies,  most  obviously 
prolate  spheroids  but  others  as  well.  Our  earlier  studies  demonstrated  that 
in  the  neighborhood  of  the  nose,  the  governing  equations  for  a  spheroid  can  be 


transformed  Into  a  form  In  which  the  paraboloid  equations  appear  as  a  natural 
limit  (t  -*■  0,  see  Eq.  (17))  in  which  the  various  parameters  of  the  flow 
remain  finite.  Hence  we  may  confidently  expect  that  any  special  feature,  such 
as  separation,  which  appears  In  the  one  will  also  appear  ,1n  the  other.  Wang 
and  Fan  (1982)  have  studied  the  unsteady  boundary  on  a  prolate  spheroid  for  ■ 
a  *  45°  (the  value  of  t  was  not  supplied).  They  found  that  whereas  the 
steady-state  boundary  layer  separated  at  p  ■  0.39/t,  the  unsteady  boundary 
layer  does  not  separate  for  any  finite  t.  This  result  we  regard  as  being  In 
direct  conflict  with  ours  and  moreover  a  matter  of  surprise.  For  not  only 
does  a  small  nonzero  t  make  very  little  difference  to  our  form  of  the 
equations,  but  the  result  Is  hard  to  understand  In  the  general  context  of 
boundary- layer  theory.  At  present  the  discrepancy  Is  Inexplicable.  The 
choice  of  coordinate  systems  for  integration  around  the  nose  does  seem  to  be 
very  cumbersome  and  It  would  be  Interesting  to  repeat  the  calculations  using 
the  system  advocated  here.  For  the  purposes  of  this  discussion  we  shall  set 
this  calculation  at  one  side. 

We  repeat,  therefore,  that  In  our  view  the  phenomena  we  have  discussed  in 
connection  with  the  paraboloid  would  also  appear  in  a  comparable  study  of 
prolate  spheroids.  In  particular  the  singularity  would  appear  after  a  finite 
time  In  an  Integration  of  the  unsteady  equations  on  the  line  of  symmetry  in 
the  neighborhood  of  the  nose  and  just  downstream  of  separation  If  o  >  ctc. 
Generalization  of  this  conclusion  to  the  unsteady  boundary  layer  over  the  whole 
spheroid  is  of  interest.  In  principle  the  integration  of  four  dimensional 
boundary- layer  equations  (x,0,y,t)  wavelike  with  respect  to  (x,0,t)  using 
the  Keller-box  scheme  Is  no  more  difficult  than  those  for  steady  two- 
dimensional  flow,  provided  we  make  the  appropriate  modifications  to  account 
for  the  varying  direction  of  flow  across  the  boundary  layer.  It  just  takes 
longer  because  there  are  more  mesh  points  to  consider  in  the  (x,0,t)  space. 
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We  may  expect  that  the  solution  of  these  equations  will  develop  a  singularity 
on  the  normals  to  a  curve  C  on  the  spheroid  after  a  finite  time*  the  precise 
time  varying  from  point  to  point  of  C,  and  It  will  be  Invlscld  In  character. 
Indeed  van  Domnelen  and  Shen  (1982b)  (see  also  Cowley  1982)  have  suggested  a 
generalization  of  their  two-dimensional  structures  which  seems  likely  to  be 
appropriate.  From  our  experience  with  two-dimensional  unsteady  boundary 
layers  we  are  confident  that  the  Keller-box  scheme  can  carry  the  Integration 
up  to  the  first  onset  of  the  singularity.  It  remains  to  be  seen  whether  the 
calculation  can  be  extended  beyond  this  time  and  over  what  part  of  the 
flowfleld  -  certainly  the  front  portion  of  the  body  can  be  treated  but  the 
precise  boundaries  are  presumably  fixed  by  an  application  of  the  Raetz  (1957) 
theory  of  Influence-regions. 

Finally,  It  Is  of  Interest  to  consider  the  possible  Impact  of  this  study 
on  the  calculation  of  steady  boundary  layers  on  bodies  of  revolution.  The 
extensive  computations  In  t  on  prolate  spheroids  showed  that  on  the  windward 
part  of  the  separation  Tine  the  flowfleld  develops  a  Goldstein-Brown  (1965) 
singularity.  On  the  leeslde  of  the  ok  of  accessibility  (Cebecl  et  al.  1981b), 
the  calculations  for  large  a  In  general  must  be  terminated  by  the  external 
streamline  through  the  ok,  so  that  the  question  of  the  structures  of  the  lee¬ 
slde  separation  line  is  Irrelevant.  At  small  values  of  a  (=6°)  this  limi¬ 
tation  does  not  apply  but  nevertheless  the  computation  breaks  down.  It  Is 
believed  that  this  breakdown  occurs  very  near  leeslde  separation  but  a  recon¬ 
ciliation  with  the  Goldstein-Brown  structure  was  not  achieved.  In  particular, 
the  blowing  velocity  is  negative  whereas  the  theory  requires  it  to  be  large 
positive. 

A  feature  of  the  singularity  Is  that  when  t  -  t  is  small  (39),  the 
blowing  velocity  Is  large  and  positive  as  x  -*■  xQ  and  after  peaking  becomes 
negative  In  order  to  bring  the  displacement  thickness  back  to  more  moderate 
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values.  Thus  there  Is  some  parallel  here  with  the  leeslde  separation.  Further 
encouragement  to  this  notion  comes  from  the  structure,  Eq.  (39),  which 
suggests  that  In  the  central  regime  of  the  singularity,  the  velocity 
components  are  almost  constant.  If  this  Is  the  case.  In  the  leeslde 
separation  singularity  the  governing  equations  for  steady  three-dimensional 
and  unsteady  two-  dimensional  boundary- layer  flow  are  effectively  the  same  and 
so  the  assumption  of  a  van  Dommelen  and  Shen  singularity  Is  as  consistent  for 
the  one  class  as  the  other. 

The  boundary  layer  on  a  paraboloid  at  Incidence  Is  of  interest  In  connec¬ 
tion  with  this  line  of  thought.  For  with  a  <  ac,  Its  structure  when  p  »  1 
has  much  In  common  with  that  of  the  unsteady  boundary  layer  on  a  circular 
cylinder,  both  on  the  leeslde  line  of  symnetry  and  on  the  windward  side.  In 
comparing  the  two  we  must  relate  p  to  t  and  the  cross-flow  velocity  on  the 
paraboloid  to  the  external  velocity  on  the  circular  cylinder  and  perhaps  we 
should  not  push  the  analogy  too  far  in  the  reversed  flow  regime  on  the  circu¬ 
lar  cylinder.  The  unsteady  boundary-layer  separates  at  e  *  106°  When  t  Is 
large  and  the  solution  Is  regular  there.  The  same  will  be  true  for  the  steady 
problem  but  it  remains  to  be  seen  whether  a  singularity  occurs  In  the  neighbor¬ 
hood  of  9  *  106°  at  finite  p.  It  is  tempting  to  conjecture  that  the  boundary 
layer  will  also  breakdown  somewhere  near  0  *  111°  with  a  singularity  at  a 
finite  value  of  p,  and  that  further  downstream  there  will  be  a  region  of  inac¬ 
cessibility,  spreading  out  In  both  directions  of  9  and  terminating  at  0  * 

106°  and  0  *  ir  as  p  +  ®,  just  as  occurs  with  the  ok  at  higher  angles  of 
attack. 
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(a)  Streamwise  skin-friction  component. 
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Figure  2.  (b)  Circumferential  skin-friction  component. 


Figure  2.  (c)  Streamwise  displacement  thickness. 


Figure  3.  Variation  of  streamwise  skin-friction  component  with  p  and  t  for  o  =  40°. 


(a)  Streamwise  skin-friction  component. 
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Figure  4.  (d)  Circumferential  displacement  thickness. 


